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SUMMARY 

A  mathematical  method  is  presented  for  the  determination  of  the  normal 
araisymmetric  end  loads  from  known  values  of  the  bulk  stress  on  a  lateral  surface  of  a  three 
dimensional  rectangular  prism.  This  extends  previous  work  which  considered  only  the 
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and  results  for  an  arbitrarily  chosen  case  are  presented. 
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SYMBOLS 


Cartesian  co-ordinate  axes 

2a, 2b, 2c 

Length,  depth,  breadth  of  the  prism 

V 

Poisson’s  ratio 

Normal  stress  components  in  the  x,y,z  directions 

2k\h,  2k2C 

Lengths  of  loaded  portion  in  the  y  and  z  directions 

^mNl  -KnI)  Llm 

Fourier  coefficients  (Loading  Antisymmetric  in  y) 

l,m,N 

Indicies  of  the  Fourier  series  (Loading  Antisymmetric  in 

y) 

Huni  ^1)  Liu 

Fourier  coefficients  (Loading  Antisymmetric  in  z) 

l,M,n 

Indicies  of  the  Fourier  series  (Loading  Antisymmetric  in 

/(y.z) 

Stress  distribution  function  on  the  faces  x  =  ±o 
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1. 


INTRODUCTION 


The  stress  distribution  in  rectangular  prisms  subjected  to  certain  normal  end  forces 
on  its  two  opposite  faces  is  a  fundamental  problem  in  Applied  Mechanics.  Reference  [1] 
supplies  a  method  for  calculating  a  symmetric  end  loading,  assuming  that  the  bulk  stress 
((Ti  +  o-y  -f  <T,)  on  a  lateral  surface  of  the  rectangular  prism  is  known  (i.e.  measurable). 
Such  a  measurement  could  be  obtained  using  advanced  thermal  emission  techniques.  The 
aim  of  this  paper  is  to  extend  this  concept  to  the  case  where  an  antisymmetric  end  loading 
is  applied. 

A  method  for  calculating  the  bulk  stresses  for  a  given  antisymmetric  end  loading  is 
presented  in  (2).  This  will  be  referred  to  as  the  Direct  solution.  In  the  Inverse  method 
the  end  loading  will  be  determined  from  a  known  distribution  of  bulk  stress  values  on  a 
lateral  surface.  In  this  case,  bulk  stresses  on  the  lateral  surface  wiU  be  simulated  using 
the  Direct  solution. 

Figure  1  shows  the  orientation  of  the  rectangular  prism  under  consideration.  The 
equations  supplied  in  [2j  allow  calculation  of  the  stress  distribution  on  a  lateral  surface 
for  an  applied  loading  which  is  antisymmetric  in  y  and  symmetric  in  x  and  z.  This  paper 
extends  this  solution  to  calculate  the  stress  distribution  on  the  same  lateral  surface  for  a 
loading  smtisymmetric  in  z  and  symmetric  in  x  and  y.  Inverse  solutions  for  both  loading 
cases  have  been  derived.  A  single  Inverse  solution  could  be  used  to  determine  a  loading 
antisymmetric  in  y  from  the  stress  distribution  on  the  face  y  =  6,  or  to  determine  a 
loading  antisymmetric  in  z  from  the  stress  distribution  on  the  face  z  =  c.  However,  two 
separate  solutions  ue  neccessary  to  determine  the  applied  antisymmetric  loading  if  only 
one  lateral  surface  is  available  (or  has  been  prepared)  for  bulk  stress  measurement  by 
thermal  emission. 


2.  DIRECT  SOLUTION 


The  bulk  stresses  on  the  surf&ce  y  =  +b  are  required  for  the  case  of  a  rectangiilar 
prism  subjected  to  certain  end  forces.  The  two  end  loading  cases  considered  are: 


( 1 )  Loading  antis3numettic  in  y,  symmetric  in  i  and  z 

(2)  Loading  smtisymmetric  in  2,  symmetric  in  x  and  y 

2.1  Loading  Antisymmetric  in  Y 
First  define  the  Galerkin  vector  P. 

F  =  iP*  +  jPv  + 

where 

Pi  =  - (“Mn*  cosh  -  {2u  +  au„a.  coth  QM„a)  sinh  x 

^  „  o:J,„  cosh  a„„a 

mrz  .  Mny 

X  cos - sm  — 

c  26 


n  ( 


P,.fe 

y32i8mh/3n,6 


f/?,iysinh^,,y  -  (2»'  + /?„,6taah/3„,6)  cosh/?„,y) 


X 


n-Kz  Inx 

X  cos - cos - 

c  a 


ilSfC 


I  M 


7f„cosh7(MC 


[ituz  cosh  7imZ  -  (2i/  +  7,j|fC  coth  fiuc)  sinh  hux]  x 


and 


firi  .  AfTTu 

X  cos - sin - 

a  26 


where 


2 

/Tiir\2 

^Afn  “ 

(it)  + 

(t) 

Pl,=  ' 

(“)  +  ( 

- 

“  / 

2 

Miry 

7,m  = 

uj  n 

26  J 

l,n  =  0, 1,2,3... 

and 

M  =  1 

(1) 
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The  stress  components  c^,  and  Ty*  are  then  related  to  F  by  the  following  equations. 


=  2(l-.)|VF,+  (.V"-£)divF, 


This  gives: 

'^*=EE 

M  n 


H  CL 

- ^ - [(1  +  ai,„acotha:jM„a)coshQM«a:  -  a*,,!  sinh  aM„a:] 

cosn  ajvrn^ 


Miry  nirz 


X  sin  • 


TItTZ  \ 

^  +  EE- 


K^,h 


2b  c  ;9*,c2sinh^„,6 

X  [( 1  -  /Jnib  tanh  /3„,5)  sinh  l3^,y  +  /3„,y  cosh  /9„ij/] 
Liuc 


2un^ir^  sinh/3, 


nirz  lirx 

cos - cos - 

c  a 


^  E  E  ^2^62  cosh  7ia,c 


M 


X  cosh  7)Af  z  +  7iiu2  sinh  ^im^]  cos  sin 

o 


2*' — - —  cosh  7ijifZ  + 
4 

Inx  .  Miry 

IT 


(2) 


(3) 


7,vCCoth7wc) 


(4) 


=  -EE 


M  n 
+  aM„xsmb  au„x 


2aifn^ccosh  auna 


(1  -  {2t/  4-  Qju„o  coth  Ou^a)]  cosh  a^nX 


Miry  .  nirz 


-sin 

25  c 


TITTZ  V  ^ 

— +EE 


n  I 


Kntnirb 
/?„,c  sinh  /3„)5 


/3„,y  sinh  I3„,y 


/3„,5  tanh  0„,b  cosh  0„,y 
7,mC  coth  ituc  sinh  7,ji,z  -  7,j^2  cosh  7ia,z 


LimMitc 


nirz  lirx  f  ^  ^ — % 

sin - cos - !•  >  /  „  .  , 

c  a  " "  27,*,5cosh7nfC 


Miry  lirx 

sin  — : —  cos - 

26  a 


(5) 


Expressions  for  <ry,<Tzj‘’'zv  '’zz  can  be  written  in  a  similar  manner. 
With  reference  to  Figure  1,  the  following  boundary  conditions  apply. 


on 

y  =  ±6 ; 

o 

II 

tT 

O 

11 

and 

r,F  = 

0, 

on 

z  =  ±c ; 

<T,  =  0, 

T*.  =  0 

and 

II 

H 

i.* 

0, 

on 

X  =  ±a ; 

< 

II 

O 

11 

and 

T«  = 

=  0. 
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AD  the  shear  stress  boundary  conditions  are  exactly  satisfied  by  (5)  and  similar  ex¬ 
pressions  for  Tiy  and  The  norma]  stress  boundary  conditions  are  approximately 
satisfied  by  equating  low  order  Fourier  terms  in  equation  (4),  and  in  similar  expressions 
for  (Ty  and  Taking  the  double  Fourier  transforms  of  these  three  equations  (with  their 
approximated  normal  stress  boundsu-y  conditions)  gives  the  expressions  (7), (8)  and  (9). 


E 

M 


J 


'2^2\ 


.  (^) 


^  ^  iLiuSi  tanh7isrc(-l)"+<^-^)/^ 
M 


7iM6(7i*fC^  +  n*ir2)2 


/  MtTcX  2  2 


=  0 


(7) 


where 


6i  =  2  when  n  =  0,  6i  =  1  when  n  ^  0 

5j  =  2  when  n  =  0  or  1  =  0,  62  =  1  when  n  /  0,  1^0 

^3  =  2  when  1  =  0,  =  1  when  1  0 


where 


tanh  aMna(-l)'+" 


asfnC(aJ,„o*  + 

2 


MnaY 


\  2b 


- 


2 


+ 


^  "  0.^c\0lP  ^ 


4  /  '  V  c  / 


-i-wi(i+-.4y  )=» 

\  8inh27,Mc/ 


(8) 


Si  =2  when  1  =  0,  =  1  when  1/0 
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Equations  (7), (8)  and  (9)  may  then  be  solved  simultaneously  for  the  unknown  Fourier 
coefficients.  The  stress  at  any  point  in  the  prism  may  then  be  calculated  by  substituting 
the  Fourier  coefficients  into  equations  (4), (5)  and  similar. 

Since  the  given  equations  are  valid  only  for  self  equilibrating  end  loads,  it  is  necessary 
to  modify  the  required  end  loading  f{y,z)  so  that  the  net  force  and  moment  acting  on 
each  end  face  is  zero.  This  produces  a  modified  loading  distribution  h(y,z). 

f2{y,z)= 

where  M  is  the  total  moment  due  to  f(y,z),  i.e. 

f{y-,^)ydydz  (11) 

The  previous  equations  calculate  the  stresses  due  to  the  modified  loading  f2{y,z).  The 
stress  component  <7*,  calculated  from  equation  (4),  is  a  Fourier  series  approximation  to 
/2(».2)- 
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2.2  Loading  Antisymmetric  in  Z 


Again  define  the  Galerkin  vector  F. 


F  =  iFr  +  jF,  +  kF, 


where 


Fi  =  V  ~z - - (qwat* cosh a„NX  -  (2i' +  a„/fa  coth  a.„Ara)  sinh  q„w*1  x 

miry  .  Nirz 

X  cos  — ; —  sin - 

b  2c 


N  I 


Ks,b 

0]tn  cosh  0Kib 


\0my  <^osh0my  -  (2i^ 


■Jf  0 mb  coth  0mb)  sinh  0my]  x 


X 


lirx 

cos - sin 

a 


Nirz 


-EL 

l  m 


^  sinh-yi^c 


7lm 


[')(i„zsinh7. 


(2i^  +  7,„ctanh7,„c)cosh7,„zjx 


and 


where 


(12) 


f,m  =  0, 1,2,3...  and  A^  =  1,3, 5,7... 
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Substituting  these  equations  into  equations  (2)  and  (3)  gives: 


m  N 


cosh  a^t^a 


[(1  +  a«//acoth  a„//a)  cosh  a„,jvar  —  ovnjvar  sinh 


X  cos 


miry  Nvz  ^ 


Ktnb 


N  I 


/3^c^cosh/3>nh 


AT^ir^  , 

2i'—^  cosh  ;JniJ/  +  1  ~  1 


]  /tTI 

X  [(1  -/3M6coth/3M^>)cosh/(?OT»  +  i9AnJ'siDl>^My]  sin cos  — 

2 

+  - f2i/m^ir^sinh')f,„2  +  (— )  ((1  - '),„ctauh7,„c) 

^  ^  smb'ri^cl  \  a  } 

1  It^x  mny 
X  sinh  7, „  2  +  7,„2  cosh7i„2|  cos - cos  —— 


(13) 


2a^^fcc  cosh 


+  q„nX  sinh 


(1  -  (2i/  +  a„jv  i  coth  a„Ara)]  cosh  cv.w* 
JVir2  .  mjry 


2c 


■  sin  - 


EE 

N  I 


2/3jviC  cosh 


/3;.,y  cosh  iIImv 


JVirz  /irx  r—v 

-  coth  Pm6  sinh  PmV  cos  — —  cos - +  >,  /  . 

Ic  a  ‘-r 

I  m 


L,„miTc 

7,„bsinh7,„c 


7,„2sinh7i„2  -  7i„ctanh7imC.  osh7i„2 


.  miry  hx 

sin  — r—  cos - 

b  a 


(14) 


Again,  expressions  for  a'y,a'i,Txy  and  t^z  can  be  written  in  a  similar  manner. 

The  boundary  conditions  of  (6)  still  apply,  with  the  exception  that  fs(y,z)  is  now  used 
to  denote  the  modified  applied  loading,  i.e. 


on 

y  =  ±b ; 

<T,  =  0, 

T^y  =  0 

and 

0 

II 

on 

2  =  ±c ; 

Oz  ~  0, 

Tj-z  ~  0 

and 

N 

II 

0 

on 

X  =  ±a  ; 

O’*  =  fs(y,  2),  Tj,  =  0 

and 

0 

11 

N 

(15) 
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Again,  all  the  shear  stress  boundary  conditions  are  automatically  satisfied  by  (14) 
and  similar  expressions  for  and  Tu-  The  normal  stress  boundary  conditions  are  ap¬ 
proximately  satisfied  by  equating  low  order  Fourier  terms  in  equation  (13),  and  in  similar 
expressions  for  Oj,  and  Oj.  Taking  the  double  Fourier  transforms  of  these  three  equa¬ 
tions  (with  their  approximated  normal  stress  boundary  conditions)  gives  the  expressions 
(16), (17)  and  (18). 


where 


4fl„.,,tanha„;„a(-l)*+*”  F  / 

^  -1  Pjr*)2  [  \  2c  / 


4ita^coth7,„c(-l)"‘+<^ 


')'lm6[7lic2  +  — 


Kt)( 


+ 


2-2  ' 


N^TT 


/  mvc 

'  4 

=  0 


6]  =  2  when  /  =  0,  6j  =  1  when  I  ^  0 


(16) 


■  tanh  a„j/<z(  — l)*"*"*^ 


(imra  'i  ^  ,  ,  ,  , 

—  j  (a^X  +  'V) 


W^),v 


N 


4Aw62tanh/?„6(-l)”>+l^->)/2 


(17) 


where 


Si  =  2  when  m  =  0,  6]  =  1  when  m  /  0 

^2  =  2  when  1  —  0,  S2  =  1  when  /  5^  0 

63  =  2  when  m  —  0  or  1  =  0,  63  =  1  when  m  ^  0,  1/0 
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a  c  ,  2a„>ra  ^Kfnt&ah 

V  sinh2a„^a  J  ^  M0l,b^  + 

’"  K^)  <A‘’ +  "•’*’)+  (t)  ’"‘*’1 


//irc\VVl  , 

*[-)  — r-'-” 


118) 


where 

^1  =  2  when  /  =  0,  =  1  when  1^0 

and 


=-r7' 

abcj^i,  J_, 


■  ma-j  .  JVjtz 

/3(3/.  *)  cos  — ^  sm  dy  dz 


Equations  (16),(17)  and  (18)  may  then  be  solved  simultaneously  for  the  unknown 
Fourier  coefficients.  The  stress  at  any  point  in  the  prism  may  then  be  calculated  by 
substituting  the  Fourier  coefficients  into  equations  (13), (14)  and  similar. 

The  required  end  loading  f(y,z)  must  again  be  modified  so  that  the  net  force  and 
moment  acting  on  each  end  face  is  zero.  This  produces  a  modified  loading  distribution 

h(y^^)- 


/3(y,z) 


/(y>2) 


ZMz 

46c^ 


where  M  is  the  totsd  moment  due  to  f{y,z),  i.e. 


(19) 


A/  =  /^  /  f{y,z}zdydz  (20) 

The  previous  equations  calculate  the  stresses  due  to  the  modified  loading  h{y,z).  The 
stress  component  o’,,  calculated  from  equation  (13),  is  a  Fourier  series  approximation  to 
fiiVjZ). 
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3. 


INVERSE  SOLUTION 


The  main  aim  of  this  section  is  to  provide  an  Inverse  solution  to  the  problem,  which 
involves  calculation  of  the  end  stress  o*  from  the  known  bulk  stress  distribution  (<Ti  + 
<Ty  +  <Ti)  on  one  lateral  surface  of  the  prism.  Consider  the  case  when  we  know  the  bulk 
stress  distribution  on  the  surface  y  =  +6.  On  this  surface,  =  0,  and  so  the  bulk  stress 
is  <Tj  +  Oj.  As  for  the  Direct  solution,  two  loading  cases  are  considered,  i.e. 

(1)  Loading  antisymmetric  in  y,  symmetric  in  x  and  z 

(2)  Loading  antisymmetric  in  z,  symmetric  in  i  and  y 


3.1  Loading  Antisymmetric  in  Y 


By  summing  equation  (4),  and  a  similar  equation  for  (t^,  the  following  expression  for 


the  bulk  stress  is  obtained. 


cosh  ciif„a 


I  ( 1  +  QAf„a  coth  an,„a)  cosh  aM„x  -  au„x  sinh  a^nX  | 


.  Miry  nvz  ^  ^  K^ib 

X  sin - r—  cos -  4-  >  >  ■  ;  - - 

2b  c  ^  ^ 


2i/n^7r^  sinh^niV  + 


|(1  -  ;^ni^tanh^„,6)sinh^„,y  +^„/J/cosh;3n,y}  cos  cos 
^  J  J  c  a 

?  ^(v) 


7i«ccoth'y,Mc) 


X  cosh  Kmx  +  1[imZ  sinh  7,j^z 


>  COS  — 

J  J  a 


Ijtx  .  Myry 

- sin  ■  — 

a  26 


co8h7,„c 


|(1  +  "rufC coth 7imc) cosh 7im;  -  7im^ sinh 7, ^zj 


Itrx  .  Miry  •r~~' 

X  cos - sin  — - - 1-  >  > 

a  7h  ^ 


a^^b^coshoMnO  [  4 


2i/ — ; —  cosh  a*,„r  + 


(=?)' 


X  I ( 1  -  aji,„a  coth  aj|/„a)  cosh  ajf,i  +  o-unX  sinh  OMn®  j  sin  cos 


+  E  E  {(1  -  ^-fctanh^„b) 

8inh/3„,y  +  /3„,y  cogh/3„,y  |  j  cos  cos  ^ 
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By  taking  the  Fourier  transform  of  this  equation,  as  shown  in  Appendix  A,  the  following 
results  are  obtained. 


abc 


(Tr  +  <rz 


Ivx  nirz 

cos - cos  — —  dx  dz 

a  2c 


r  ,  If  Af2jr2 

2 

aL  + wl*"— 7— 

^(=?y(7y}i- 


r  3  ,  If 

2  /n7r\2' 

J 

+ 


+  cH'^)  (2;/  +  1 


-  /3„,6tanh  /3„|5  +  /?„,6coth  /3„|6^ 


(22) 


where 


=  2 

when 

n  =  0, 

^1  =  1 

when 

n  0 

=  2 

when 

n  =  0 

or  I  = 

0,  62 

=  1  when  n  ^  0,  1  7^  0 

=  2 

when 

1  =  0, 

^3  =  1 

when 

Equation  (22)  may  then  be  applied  for  each  value  of  n  and  I  to  give  a  set  of  linear 
simultaneous  equations  with  the  Fourier  coefficients  K„i  and  L/m  as  the  unknowns. 
Since  the  boundary  conditions  of  equation  (6)  stiU  apply,  equations  (7)  and  (8)  are  also 
used  in  conjunction  with  (22).  In  effect,  equation  (22)  replaces  equation  (9)  for  the  Inverse 
solution.  This  imposes  the  restriction  that  1  =  (Af  -  l)/2,  so  that  the  number  of  linear 
equations  produced  by  (7)  (8)  and  (22)  equals  the  number  of  unknown  Fourier  coefficients. 
The  resulting  set  of  equations  is  then  solved  to  produce  a  set  of  Fourier  coefficients  Htf„, 
Kni  and  Lm  for  the  Inverse  problem.  Equatioii  (4)  is  then  used  to  obtain  the  applied 
stress  ffi- 


5 

f' 


n 


3.2  Loading  Antisymmetric  in  Z 


By  summing  equation  (13),  and  a  similar  equation  for  the  following  expression  for 
the  bulk  stress  is  obtained. 


<^6  =  ^  a  { ^  ^  “mwu)  cosh  a,„jsri  -  a„f,x  sinh  | 


^  cosha„N<i  1 
mxp  .  Nirz 


miry  .  jynz  V''  Ajno  ivi 

X  cos  — p-  sin  — - (-  >  >  -z-z - p-— ,  21^ - 

*  2c  ^^c^cosb^^bl  4 


iiTmfc  L  •'VV  ,  ^  /iTTc^ 

__|2.__eosh^^y+(-j 


X  I  ( 1  -  coth  b)  cosh  ^mV  +  ^mV  sink  j  |  sin  ^  cos  ^ 

+  EE  ^ +  (— )  {(1  -  7,.ctanh-n„c) 

.1  ill  ix’*  mir« 

X  sinh  7,„  r  4-  7,™  2  cosh  7i„  2  >  cos - cos  -7— 

i  j  a  b 

+  E  E  { ( 1  +  7to  c  tanh  7,,,,  c)  sinh  7,„,  z  -  7,„  2  cosh  7,„  z  j 

litx  mzy  flmN®  f„  2  2  ,  ( Nzb\‘ 

X  cos  cos  +  z.  2..  'i  a - r -  2uTn‘z^ cosh a„Ar*  +  I  -;p'  ) 

®  “‘^«mw*^cosha,„;vol  \  2c  7 

X  I  ( 1  -  amN<i  coth  a„,K<i )  cosh  a„,yvx  +  q;„,/,x  sinh  a„;,x  1 1  cos  sin 
+  E E -^/~^0-b  [2-'^-^cosh/3„y  4-  {(1  -0„bcoth0„b) 


X  cosh  0my  +  0Nty  sinh  /3jvjy  |  sin  cos 


By  taking  the  Fourier  transform  of  this  equation,  as  shown  in  Appendix  B,  the  following 
results  are  obtained. 


Iwx  .  Nirz 

<Tt)\  .  COS - Sin  — - —  ax  az 

a  2c 


=  E 


+E 


2  1““''  '  a„;,62 


m 


tanh  a„Na 


6[7,i  +  (^f] 


(“)’(^)1] 


coth7i„c 


^2,a2c2 


+  ^21/ +  1 — jSji/jb coth /3f/ib -h (24) 


where 

=  2  when  1  =  0,  =  1  when  /  0 

As  for  the  case  of  an  applied  loading  antisymmetric  in  3/,  a  set  of  Fourier  coefficients 
ftnw,  Km  and  Li„  may  be  found  by  solving  the  set  of  simultaneous  equations  produced 
by  (16),(17)  and  (24).  Equation  (13)  is  then  used  to  obtain  the  applied  stress  (t*. 
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4. 


RESULTS  AND  DISCUSSION 


The  Direct  and  Inverse  programs  were  implemented  using  FORTRAN  77  on  an  Apollo 
DNlOOOO  computet.  The  Gaussian  Quadrature  method  was  used  for  the  integrations 
required  by  equations  (9),(18),(22)  and  (24).  The  resulting  sets  of  linear  simultaneous 
equations  were  solved  using  Gaussian  elimination  with  full  pivoting.  Double  precision 
was  used  for  all  calculations. 


4.1  Direct  Results 


The  equations  given  in  the  previous  sections  are  applicable  to  a  prism  with  an  arbitrary 
rectangular  cross  section.  However,  for  illustrative  purposes,  a  square  cross  section  with 
a  square  loading  area  was  considered,  where  b  —  c  =  2  and  ki  =  k2  —  0.5.  First  consider 
a  loading  symmetric  in  x  and  y,  and  antisymmetric  in  z.  With  reference  to  Figure  3,  a 
stress  of  f{y,z)  =  ±16  was  applied  inside  the  loading  area,  with  /(j(,r)  =  0  outside  the 
loaded  area.  Equations  (19)  and  (20)  must  first  be  used  to  produce  a  self-equilibrating 
loading  f3(y,z)  on  the  end  face.  FVom  (20)  the  total  moment  is: 

/+C 

J  f{y,z)zdydz 


t>  +  k7C 

/  /  —l^zdz  + 

1  \6zdz 

1 

1 

Jo 

dy  +  0 


=  16fc2C*2fcj6 
=  32klkibc'‘ 
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Applying  (19)  gives  the  modified  loading  function. 


Hence, 


/3(k,z) 


/(l/.z) 


=  /(y.z)- 


ZMz 

4ic^ 

24fc|{;iz 


=  /(y,z)  -  l-5z 


h(yj^}  —  —1-52  outside  the  loading  area 

=  +16  —  1.5z  inside  the  loading  area,  0  <  z  <  k2C 
=  —16  —  1.5z  inside  the  loading  area,  —  k2C  <  z  <  0 

Sample  residts  from  the  Direct  method  for  an  applied  loading  antisymmetric  in  z  are 
shown  in  Figure  4.  Plots  of  the  bulk  stress  distribution  on  the  y  =  +b  surface  and  the 
end  stress  o*  are  presented.  These  results  were  obtained  using  seven  terms  of  the  Fourier 
series  (i.e.  /,m  =  0, 1,  ...7  and  N  =  1,3,  ...15),  which  showed  good  convergence.  Shown  in 
Figure  4a  is  the  surface  plot  of  o*,  with  its  corresponding  values  indicated  in  the  contour 
plot.  Figure  4b.  The  plots  show  the  stress  values  ranging  from  -18.0  to  +18.0,  in  good 
agreement  with  the  applied  stress  field  /j(y,z).  The  bulk  stress  distribution  cri,  on  the 
3/  =  +6  surface  is  shown  in  Figures  4c  and  4d. 


Similar  results  are  shown  in  Figure  5  for  an  applied  loading  which  is  antisymmetric 
in  y,  but  symmetric  in  x  and  z.  The  modified  loading  function  /2(j/,z),  produced  by 
equations  (10)  and  (11)  is  shown  in  Figure  2  as: 

hiVi^)  =  ~l-5j;  outside  the  loading  area 

=  +16  —  \.by  inside  the  loading  area,  0  <  y  <  kjb 
=  —16  -  1.5y  inside  the  loading  area,  —  kjb  <  y  <  0 


16 


4.2 


Inverse  Results 


In  the  Inverse  problem,  the  bulk  stress  field,  generated  by  the  Direct  method,  is  used 
to  calculate  the  applied  end  stress.  For  comparison,  the  results  of  the  stress  distribution 
(Tx  obtained  from  both  the  Direct  and  Inverse  methods,  at  different  Fourier  indicies,  are 
presented.  Figure  6  shows  the  results  obtained  for  an  applied  loading  antisymmetric  in  z, 
symmetric  in  x  and  y.  Figure  7  shows  similar  results  for  an  applied  loading  antisymmetric 
in  y,  symmetric  in  x  and  z. 

Reference  [Ij  presented  similar  results  for  the  stress  distribution  for  a  prism  with 
symmetric  end  loads.  It  also  demonstrated  that  more  Gauss  points  are  needed  for  numeric 
integration  to  obtain  a  reasonable  Inverse  solution  for  a  high  number  of  Fourier  coefficients. 
The  8mtis3rmmetric  solution  showed  a  similar  trend.  The  results  presented  in  Figures  6 
and  7  use  the  minimum  number  of  Gauss  points  to  obtain  a  reasonable  solution.  Further 
increases  in  the  number  of  Fourier  coefficients  and  Gauss  points  were  not  warranted  for 
preliminary  work  due  to  the  large  increase  in  computer  resources  required. 


5.  CONCLUSION 

This  work  has  successfully  extended  previous  work  on  the  stress  distribution  in  rect¬ 
angular  prisms  to  include  the  case  of  antisymmetric  end  loading.  This  distribution  may 
be  calculated  from  a  set  of  bulk  stress  measurements  on  one  lateral  surface  of  the  prism. 
Analysis  has  been  performed  for  the  case  where  the  bulk  stress  distribution  used  in  the 
Inverse  method  was  produced  by  the  Direct  method. 

Although  the  concept  has  been  successfully  proven,  a  considerable  amount  of  work 
will  be  required  to  extend  this  methodology  to  a  procedure  which  can  be  routinely  used 
for  the  analysis  of  thermoelastic  measurements  from  real  structures. 
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FIG.  1  RECTANGULAR  PRISM  WITH  CENTRAL  LOADING 


FIG.  2  APPLIED  LOADING  ANTISYMMETRIC  IN  Y 


<y<  k,h 
“  <  y  <  0 


<  «  <  fcjc 
-  Arjc  <  ;  <  0 


FIG,  3  APPLIED  LOADING  ANTISYMMETRIC  IN  Z 


Fig.  4a  Applied  Eod  Loading 


Fig.  4c  Bulk  Streia  on  y 


Fig.  4b  Applied  End  Loading  Fig.  4d  Bulk  Stieu  on  y 


FIG.  4  DIRECT  SOLUTION  FOR  APPLIED  LOADING 
ANTISYMMETRIC  IN  Z 


Fig.  5a  Applied  End  Loading 


Fig.  5c  Bulk  Stress  on  y  =  6 


Fig.  6b  Applied  End  Loading 


Fig.  6d  Bulk  Stress  on  y  = 


FIG.  5  DIRECT  SOLUTION  FOR  APPLIED  LOADING 
ANTISYMMETRIC  IN  Y 


Direct  Solution 


Inverse  Solution 


Fig  6a.  l,m  =  0,1;  =  1,3.  Two  point  integration. 


Direct  Solution  Inverse  Solution 

Fig  6b.  l,m  =  0, 1,2;  iV  =  1, 3, 5.  Two  point  integration. 


Direct  Solution  Inverse  Solution 

Fig  6c.  I,m  =  0, 1...3;  JV  =  1,3. ..7.  Three  point  integration. 


Direct  Solution 


Inverse  Solution 


Fig  6d.  1,  m  =  0, 1...4;  JV  =  1,3...9.  Four  point  integration. 


Diiect  Solution  Inverse  Solution 

Fig  6e.  l,m  =  0,1...5;  N  =  Four  point  integration. 


Direct  Solution  Inverse  Solution 

Fig  6f.  I,m  =  0, 1...6;  N  =  1, 3. ..13.  Five  point  integration. 


Direct  Solution  Inverse  Solution 

Fig  6g.  I,m  =  0,1...7;  JV  =  1,3. ..16.  Five  point  integration. 


FIG.  6  DIRECT  AND  INVERSE  SOLUTIONS  FOR  APPLIED 
LOADING  ANTISYMMETRIC  IN  Z 


Du«ct  Solution  lovene  Solution 

Fig  7u.  {,n  =  0, 1;  Jif  =  1,3.  Two  point  integtntion. 


Dit«ct  Solution  Inveise  Solution 

Fig  7b.  {,n  =  0, 1,2;  Af  =  1,3, 6.  Two  point  integtation. 


Direct  Solution  Inverse  S<Jntion 

Fig  7c.  (,n  =  0,1. ..3;  M  =  1,3.. .7.  Three  point  integration. 


Direct  Solution  Inverse  Solution 


Fig  7d.  (,n  =  0, 1...4;  M  =  l,3...l).  Four  point  integration. 


Diiect  Solntion  Invetie  Solution 

Fig  7e.  i,n  =  0, M  =  Four  point  integration. 


Direct  Solution  Inverse  Solution 

Fig  7f.  (,n  3  0, 1...6;  M  =  1,3.. .13.  Five  point  integration. 


Direct  Solntion  Inverse  Solution 

Fig  7g.  (,n  =  0,1. ..7;  M  =  1,3.. .16.  Five  point  integration. 


% 

FIG.  7  DIRECT  AND  INVERSE  SOLUTIONS  FOR  APPLIED 
LOADING  ANTISYMMETRIC  IN  Y 


X 


DERIVATION  OF  EQUATION  (22) 


APPENDIX  A 


From  Equation  (21) 


—  - 1(1  +  coth  a^^a)  coshaif^z  —  sinhQ*,„x| 


Miry  nirz 


2im*T^  sinh/%, 


iw  Ty  niTz  ^  %  K^iv 

X  sin  — 7-^  cos - +  >  >  ,  ^  - 

2b  c  0;,c^^oh/3„,b 

X  {(1  -fl.ifctanhflj6)sinh/3.,y  +  /l„vcosh/l.,!/||  cos  cos  ^ 

frt  M^x^  / lxb\^  ( 

+  ^  L  +  i  V  j  {(1-Wcothw) 

X  cosh  7,„2  T ,m2  sinh  7,^2 1 1  cos  sin 

a  2d 


1  M 

Ixx  .  Mxy 


Hu.a 


cosh  Oms^o 


IXX  M  Xy  r  ^  *-  A  i 

X  cos - sin  — ~  +  >  >  T-TT 

a  2b  ^  ^  oi.fr’ 

Af  n  *•" 

X  {( 1  -  Q^„a coth Ok.o) cosh<»j,„x  +  a^„x sinh  o„„x} 


M^x^ 

2i^ — - —  cosh 


/nir6\^ 
,..x+(— ) 


Mxy  mtt 

sin  — .  cos  - — 
26  c 


k:.,6 


X  sinh/3„,»  +  fl.tycoshj3,,y|j 


2i/l’»’5inh0„,y  +  (^)  {(1  -  Aiktanh/?,,!)) 


nx2  /»! 

cos - cos - 

c  o 


On  a  lateral  surface  of  y  =  6,  iTj,  =  0  and  so  the  bulk  stress  on  this  surface  becomes 
(Ti  +  Oj.  Putting  y  =  6  into  the  above  equation,  and  noting  that  sin  ^  =  (  — 1)**^ 
for  M  =  1,3,5...,  yields  the  following  result. 


<^b  = 
-i)/2 


APPENDIX  A 


y  - — r - >  ( 1  +  aMwOcoth  OMnO)  cosh  a,«„i  —  a**..*  sinh  a 

^Tf  cosh  cka^.a  l 

M  n  •" 

+  EE^ij2f;h^_,4  j2.m^ir^nnhA.t+  {(1  -  ^.,^tanh/3,.,^) 

.  .  -  .  .  -  "1 1  nrz  hex 


}nirz 
cos - 


X  sinh;3.,i  +  /iiicosh/SLiij 

nir2 

cos - - 

c 

,  v-v- W(-i)<"-'>^’ 

2i* 

1  M  'l'iM*’cosh->,MC 

4 

X  cosh 7, m2  +  7iM2smh7,M2 ll  cos 

/J  a 

+  ^  5  +  y<uccoihy. 

V- L  W’*-’ 

+  5  ?  f  "-4- 


nit^  ,  f, 

-^coshwz  +  (  —  j  |(1-7,mCC 

lirz 


V- L  W’*-’ 

«i>’coslia^,<i  p  4 
<  cosh  Omh^  +  <3ci4n^  sinh  «*,«* 1 1  cos 

'i 

t"  52  52  Vj  3^  ' L  [21'/’ ir’  sinh /3,, 6 

“  V /??,«’«<>»>  A.*  [ 

■»  1  n»T 


,»,c)  cosh  i,mJ  -  sinh  1,^2  )  cos  — 

a 


r’lr’  /n»fc\*r 

—  cosha«.x+  ^  — j  |(1 


-  ai,.acothau„a) 


<  sinh  (3, .16  +  /3,., 6  cosh  j  cos  -~  cos 


To  determine  the  Fourier  coefficients  Hunt  A„i  and  ijjvf  it  is  necessary  to  take  the  double 
Fourier  transform  of  the  above  equation. 


/:/:<— 


pKz  <iiez  ,  ^ 

cos - cos  - —  dx  dz 

a  c 

f+<  /+«/• 


=  f  J  [(I)  +  (U)  +  (in)  +  (lV)  +  (V)  +  (Vl)  cos^cosi—dxdc 


Note  that  p  and  q  are  now  fixed  positive  integers. 
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Integrating  the  RHS,  tahing  term  (I) 

I  I  (I)  cos - cos - dxdz 

J  —  C  J  — «  ^  ^ 

L +.^.,.otha..a)  cosh.., x  cos  <i. 
coshaji,,a  ./-a  ® 

/■+“  ,  p»* ,  1  /■■^‘  J 

-  a.,  /  *  Mnh  a.-*  <=os - “®  /  - “'• 

“  y_,  o  jy_,  c  c 


=E 


coshaM^a 


(1  +  a^nO  coth  or*(„a) - ^  sinhoi^nQ 


2a2„“(-l)'’  u 

- ! - —=  cosh  a. 


.a  + 


2ojk„( 

-!)>’[< 

— 

*ttn 

(?)’] 

ai,  + 

(?)' 

]■ 

sinh  ajf.,a 


^  2ff.„av„6iac(-l)»‘'K*'-*>^^ 


( 1  +  aMnO  coth  aMnO)  tanh  0[M..a 


(?) 


ianhaM„a 


Put  I  =  p  and  simplify  to  give  : 

V-  4ff..,af,„oc«i(-l)'+<"-‘»’  ,  ^ 

=  >  — — — - h; - 1__ - tanh  o..,o 


where 


6i  =  2  when  n  =  0,  =  1  when  n  /  0 

Taking  term  (II) 

/■+'  /■+•  P*I  qirc  ,  . 

/  /  (Ii)cos - cos - OlOi 

/-c  2—  »  « 
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nir«  9ir5  , 

cos - cos  - —  dz 

c  c 


"y32,c«nh^, 
iCioMs 


0lc 


{(1  -  /Sl..htaiih  A,J.)  +  A,hcoth  A-t} 
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where 

82  =  2  when  n  =  0  or  i  =  0,  62  =  1  when  n  ^  0,  1^0 


Taking  term  (III) 


fym)ccs?^cosSIldzdz 
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X  /  cos - cos - dx 

y-a 
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(2t/ — - 

-5  cosh  i.vc 
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V  a  /  V  a  ;  7?„+(?)’ 

Put  71  =  g,  re-arrange  and  simplify  to  give; 


■^4£,„<3ac(-l)"^-<*^-»>^^f  (lKb\ 


=E 


\-7 


(¥)' 


tanh  7i«c 


where 


^3  =  2  when  1  =  0,  =  1  when  1  /  0 
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The  integrals  of  parts  (IV)  (V)  and  (VI)  were  derived  using  the  same  approach  as  shown 
in  (I)  (II)  and  (III). 
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where 

=  2  when  n  =  0,  6i  =  1  when  n  /  0 

^2  =  2  when  n  =  0  or  1  =  0,  62  =  1  when  n  ^  0.  /  7^  0 

^3  =  2  when  /  =  0,  ^3  =  1  when  I  ^  0 
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The  individual  integral  components  may  now  be  summed  to  give  the  following. 
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Simplifying  this  equation  gives  the  following  result,  as  presented  in  equation  (22). 
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DERIVATION  OF  EQUATION  (24) 


From  Equation  (23) 
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On  a  lateral  surface  of  y  =  6,  <Ty  =  0  and  so  the  bulk  stress  on  this  surface  becomes 
crj,  =  <Ti  +  cTj.  Putting  y  =  b  into  the  above  equation,  and  noting  that  rosmir  =  (  —  I)"* 
for  m  =  0,1,2...,  yields  the  following  result. 
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To  determine  the  Fourier  coefficients  HmNi  fiNi  and  Lim  it  is  necessary  to  take  the  double 
Fourier  transform  of  the  above  equation. 
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Note  that  p  and  Q  are  now  fixed  positive  integers. 
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Integrating  the  RHS,  taking  term  (I) 
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Put  I  =  p  and  simplify  to  give  : 
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Taking  term  (III) 


-W  I(2i/m^ir^  +  (  —  ]  (1 -ii„ctanh-rtac)} 

-  2^  2^  sinh  7(me  [  I  V  “  / 

/•+'  .  <?>ri  j  /I’rfc'l’  /■"'. 

X  /  sinh^fim^sin a2  +  1  —  j  ^im  J 


cosh  Urn 


2c 


l»I  P»I 

X  /  cos - cos - ai 

J-a  “  “ 


_  EimM£(_jr.  I  f  2t,„ +  (  — )  (1  -  7rmCtanh7i„c)} 
7L*>’ '‘"h^lmC  1  ^  Vo/ 


rf-MS)  ’“*(^1 


=i; 


-(v) 

2Li„i,ac(-ir^ 


Tim 


+  (^)  ] 

X  coth7lmC 

Put  N  =  Q,  re-arrange  and  simplify  to  give; 
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where 


6,  =  2  when  /  ^  0,  6^  =  1  when  I  f  0 
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The  integrals  of  parts  (IV)  (V)  and  (VI)  were  derived  using  the  same  approach  as  shown 
in  (I)  (II)  and  (III). 
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The  individual  integral  components  may  now  be  summed  to  give  the  following. 
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Simplifying  this  equation  gives  the  following  result,  as  presented  in  equation  (24). 
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